Meissner effect and reentrance phenomenon in heated rotating nuclei 



o 
o 

(N 

c 



D.J. Dean, 1 K. Langanke, 2,3 H.A. Nam, 4 and W. Nazarewicz 5 ' 6 

1 Physics Division, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, TN 37831 
GSI Helmholtzzentrum fiir Schwerionenforschung, Planckstr. 1, D-64291 Darmstadt, Germany 
1 Institut fiir Kernphysik, Technische Universitdt Darmstadt, 
Schlossgartenstr. 9, D-64291 Darmstadt, Germany 
'* National Center for Computational Sciences, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, TN 37831 
Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996 
^Institute of Theoretical Physics, Warsaw University, ul. Hoza 69, 00-681 Warsaw, Poland 

(Dated: January 6, 2009) 

Rotational motion of heated 72 Ge is studied within the Shell Model Monte Carlo approach. We in- 
vestigate the nuclear pairing correlations associated with ,P Cooper pairs as a function of the angular 
momentum polarization and temperature. The reentrance of pairing correlations with temperature, 
a specific feature of the final system, is predicted at high rotational frequencies. It manifests itself 
through the anomalous specific-heat behavior. 

PACS numbers: 21.10.-k, 21.60.Cs, 21.60.Ka, 27.50. +e 
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The Meissner effect is a dramatic demonstration of the 
superconductivity phenomenon. When a metal is cooled 
below its superconducting transition temperature T c in 
the presence of an external magnetic field, the magnetic 
flux inside vanishes For type I superconductors, mag- 
netic flux is expelled below a critical field H c , while at 
H=H C , a sharp transition to a normal conductor takes 
place. In type II superconductors, this transition is grad- 
ual, and it is characterized by two critical temperatures, 
T cl and T c2 . For T cl < T < T c2 , the system exists in the 
vortex (mixed) state in which the normal and supercon- 
ducting phases coexist. At a fixed temperature, the flux 
penetration exists between critical magnetic fields H c \ 
and H C 2- The critical fields depend on temperature and 
pressure, and in some cases superconductivity survives in 
extremely high values of H 
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Atomic nuclei exhibit a variety of emergent phenom- 
ena, including that of fermionic pairing at low en- 
ergies, the wave functions of an open-shell even-even nu- 
cleus can be viewed as a condensate of nucleonic Cooper 
pairs. Pairing plays a crucial role in our understanding 
of low-energy nuclear structure: (f ) it is responsible for 
the energy gap in the spectra of even-even systems; (2) 
it gives rise to the odd-even staggering of nuclear bind- 
ing energies; and (3) it reduces the nuclear moments of 
inertia from the rigid body value [ESQ- The isovec- 
tor J=0 pairing, i.e., the coupling of like nucleons to 
J=0 pairs, has been identified as the dominant pairing 
mode. With increasing temperature, nuclei undergo a 
distinct transition from the pairing phase to a normal 
phase. It is expected that around &;T=0.5-f MeV [1] the 
static condensate vanishes. However, as nuclei are finite 
systems, this transition is smeared out by thermal and 
quantum fluctuations. Most theoretical studies of the 
superconducting-to-normal transition have focused prin- 
cipally on the relationship between pairing correlations 
and an associated peak in the specific heat C v — dE / dT 



(see, e.g., 0, [l(| EH Experimental investigations 

of the level densities also indicate the possibility of the 
thermal melting of the nuclear pairing phase [13J. 

The similarity between the Meissner effect in super- 
conductors and the pairing response to fast nuclear rota- 
tion was noticed in the f 960s by Mottelson and Valatin 
[Til ]. Based on the analogy between the magnetic field 
and rotational frequency cj, they predicted a sudden 
pairing collapse at high angular momenta (the so-called 
Coriolis-antipairing effect). However, it was realized a 
decade later |15l.[l6j] that the competition between the nu- 
clear pairing and rotation is more complex: the Coriolis 
force acting on nucleons occupying high-j orbits gradu- 
ally breaks Cooper pairs, thus giving rise to an increased 
rotational alignment (angular momentum polarization) 
and decreased pairing. That is, nuclei are in fact type II 
superconductors in which the vortex state can be associ- 
ated with a many-quasiparticle configuration containing 
broken nucleonic pairs [l5j |. 

Recently, there has been increased interest in proper- 
ties of polarized (asymmetric) Fermi systems with un- 
usual pairing configurations. This includes the exis- 
tence of superconductivity in a ferromagnetically ordered 
phase [13, 0, [H, H3| and an interplay between pairing, 
spin polarization, and temperature in condensates having 
an unequal number of spin-up and spin-down fermions 
[2ll [22I I23I ]. While most of the underlying theoretical 
discussion has been based on schematic models [l9|, 
or mean-field approaches (BCS or extended density func- 
tional theory) [2ll [H| , there have been very few studies 
for finite systems (nuclei, atomic condensates, metallic 
grains), based on realistic interactions in large configu- 
ration spaces, that properly take into account quantum 
fluctuations. In the context of atomic nuclei, the major- 
ity of such investigations have been restricted to rotating 
cold nuclei or to non-rotating heated nuclei, thus explor- 
ing only limits of the (u),T) phase diagram. In Ref. [25| . 
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the Shell Model Monte Carlo (SMMC) approach was 
used to study the spin distribution of nuclear levels and 
ground-state (g.s.) thermal moments of inertia of nuclei 
from the A~60 mass region. In our previous paper [Til ] . 
we explored the competition between shape deformations 
and pairing correlations as a function of temperature in 
7V=40 isotones. The generic features of pairing corre- 
lations in sd-shell nuclei as a function of excitation en- 
ergy, angular momentum, and isospin were investigated 
in Ref. [26[ within the nuclear shell model. 

The objective of this Letter is to explore the com- 
petition between symmetry-breaking modes in finite 
fermionic systems. To this end, we study the nuclear 
behavior as a function of temperature and rotation. We 
perform detailed calculations for a medium-mass nucleus 
72 Ge that is known to exhibit g.s. pairing correla- 
tions and significant quadrupole collectivity. 

We utilize the SMMC method [27[ that allows stud- 
ies of nuclei at finite temperatures with the relevant de- 
grees of freedom included; the method properly accounts 
for the thermal and quantal fluctuations that are impor- 
tant to describe phase transformations in finite-size sys- 
tems. Within SMMC, a finite-temperature observable is 
given by a thermal average. For certain classes of residual 
nucleon-nucleon interactions [28|, such as the force em- 
ployed in this work, the evaluation of observables is ex- 
act, subject only to statistical errors related to the Monte 
Carlo integration. The details of our SMMC calculations 
follow those of Ref. [ll|. Namely, the single-particle (s.p.) 
Fock space corresponds to the complete (0/lp — 0gld2s) 
shells for protons and neutrons. The nucleus 72 Ge is 
described by 12 valence protons and 20 valence neutrons 
outside the closed core of Ca. As we are concerned here 
with a description of collective quadrupole and pairing 
correlations at relatively low energies, we have employed 
a pairing+quadrupole-quadrupole Hamiltonian that was 
chosen to reproduce the low-energy spectrum of 64 Ni and 
64 Ge. In order to generate the angular momentum po- 
larization, we consider the routhian = H — ljJ z , 
where the cranking frequency uj (in units of MeV) en- 
ters through the cranking term. (For early applications 
of SMMC to high-spin states, we refer the reader to 
Refs. [Hill.) 

The thermal pro perties of 72 Ge were discussed in the 
previous paper [ll(. As one of the striking results, the 
SMMC calculations gave clear evidence for the breaking 
of isovector pairs at temperatures around kT c « 0.6 MeV 
in 72 Ge as reflected by a noticeable peak in the specific 
heat C v . (The deformed-to-spherical transition results 
in a broad shoulder in C v only.) In a model of non- 
interacting particles, the neutrons in 72 Ge would com- 
pletely occupy the fp shell. However, correlations in- 
duced by the residual interaction make it energetically 
favorable to scatter neutrons across the 2.5 MeV iV=40 
shell gap. Figure [T] shows single-neutron occupations in 
the many-body wave function of 72 Ge as a function of 
rotational frequency at two temperatures: T=0.47MeV 
(slightly above g.s.) and 1.6 MeV (well above T c ). In the 



g.s. configuration, the total neutron occupation of the 
gds shell is about 3.5, with about 3 neutrons in the <?g/ 2 
orbital. As protons only occupy about half of the fp 
shell, their excitation into the gds shell is significantly 
smaller (only about 0.3 protons are promoted). Upon 
rotating the nucleus, the s.p. cranking term u>j z , repre- 
scnting the combined effect of the Coriolis and centrifugal 
force [31j . generates angular momentum polarization by 
lifting the magnetic m-degeneracy of s.p. states. This 
is clearly seen in Fig. [lja), especially for the <?g/2, which 
has the largest j-value in our s.p. space. For instance, 
at w=0.5MeV, the occupation of m=9/2 orbit grows to 
~0.8, while the one of m —-9/2 is reduced to ^0.15 in 
comparison to the value of 0.3 at lj—Q. 
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FIG. 1: (Color online) Neutron occupations of the individ- 
ual fr/2, fs/2> an d <?g/2 S -P- orbits as a function of cranking 
frequency u> in the wave function of 72 Ge for temperatures 
fcT=0.47MeV (a) and 1.6 MeV (b). The angular momentum 
alignment in the nucleus is generated by the gradual increase 
of occupations of positive-m orbits (solid lines) at the expense 
of negative-m states (dashed lines). The doubled magnetic 
quantum numbers, 2m, of individual orbits are marked. 

The highly asymmetric pattern of s.p. occupations seen 
in Fig. UK a) is due to an interplay between the Coriolis 
force that tries to align individual single-particle angu- 
lar momenta along the axis of rotation, thus introduc- 
ing spatial polarization of the system and the symmetry- 
restoring pairing force. Indeed, in the J—0 pair operator 

At = E J ; m >o(- i r m 4m a ]-m' a11 magnetic substates 
appear with the same weights. Thus one observes in 
Fig. Q] that, at low temperatures, the occupations of the 
various m substates of the gg/% and f§/2 orbitals do not 
follow a simple thermal ordering. In fact, the occupations 
of the substates with largest m values are largest, as ex- 
pected from energy considerations, but far from one (the 
extreme s.p. limit). On the other hand, the states with 
the lowest negative m values have very similar reduced 
(but nonzero) occupations, indicating that the residual 
population, although disfavored by the cranking term, 
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allows the system to gain energy through pairing. The 
quadrupole interaction, although it also gives rise to some 
rearrangement of s.p. occupations, is less important in 
this discussion. Indeed, the total occupation of the neu- 
tron d 5 / 2 shell is about 0.4, thus indicating that the low- 
est deformed 7r=+ Nilsson orbits have a fairly pure 179/2 
character. 

If the nucleus is heated above T c , or if the rotation 
is rapid, pairing correlations are dramatically reduced 
(see Fig. [5]). However, as expected for a finite system, 
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FIG. 2: (Color online) Neutron (top) and proton (bottom) 
J=0 pairing strength (A^A) as a function of temperature 
for 72 Ge at four values of lo (in MeV). At higher rotational 
frequencies, neutron pairing locally increases at low temper- 
atures reaching a maximum (marked by arrows). 

the superconducting-to- normal transition does not occur 
abruptly but is smeared out [lj, [32[ . 

Since at large temperatures pairing and quadrupole 
correlations are quenched, the occupations of s.p. states 
should follow the Thomas-Fermi ordering governed by 
s.p. routhians, e" m m ej — Lum. This is precisely what is 
seen in Fig. [TJb). 

The expectation values of angular momentum compo- 
nents as a function of temperature are shown in Fig. [3] 
At low values of u, the angular momentum alignment 
(J z ) = (J z ) is around 2.5 Ti at the lowest tempera- 
tures considered, i.e., the main contribution to the ther- 
mal average comes from the collective superconducting 
g.s. band. At higher temperatures, the alignment in- 
creases to ~4.5 7i, reflecting the large contribution from 
the rotationally aligned {gg/2) 2 configuration involving 
one broken neutron pair coupled to J=8. At w=0.5MeV 
and low T (Fig. Wp), the lowest routhian is given by 
the aligned (39/2) 2 configuration carrying the total align- 
ment of ~10.5 Ti (=2.5+8). As the nucleus is heated 
up, the contributions from less-aligned states become 
more important and (J z ) is gradually reduced. The 
average angular momentum (J), defined through the 
relation (J 2 ) = (J)((J) + 1), steadily grows with T 
for low and high values of w. At low rotational fre- 
quencies, this increase comes from both the parallel, 

(^z)r.m.s. = y (Jz), and perpendicular angular momen- 



tum (J_i) r . m .s. = y {J 2 — J 2 )- The situation is differ- 
ent at w=0.5MeV. Here, (J z ) r .ni.s.~10-5/i and the tem- 
perature dependence of (J) comes primarily from quan- 
tum and thermal fluctuations in the perpendicular di- 
rection. The average moment of inertia, defined as 
(J) / \d{E) I 'd(J)) , steadily grows with to at low tempera- 
tures, in accordance with weakening pairing correlations. 
At fcT=1.6MeV, the moment of inertia is fairly constant 
as the static pairing is gone. 
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FIG. 3: (Color online) Angular momentum in 72 Ge as a func- 
tion of kT at cj=0.25 MeV (a) and 0.50 MeV (b). Plotted are: 
average angular momentum (J) (solid line); alignment (J z ) 
(dashed line); (Jx) r . m . s . (dash-dotted line); and (J z ) (dotted 
line). 

An interesting feature seen in Fig. [2] is the local in- 
crease of neutron pairing with temperature at low val- 
ues of T and high rotational frequencies. The growth 
of the pairing correlations with temperature was pre- 
dicted in Ref. [33| and discussed in Ref. [l9j ]. based on 
a schematic seniority model. This phenomenon is a con- 
sequence of the quantum fluctuations in a finite system. 
At T«0 and large values of ui, the yrast line corresponds 
to a strongly aligned, many-quasiparticle configuration 
with strongly reduced pairing. With increasing temper- 
ature, less-aligned excited configurations with lower se- 
niorities enter the canonical ensemble, and this reintro- 
duces the pair correlations. At still higher temperatures, 
the Meissner transition takes place and pairing correla- 
tions decrease. The reappearance of pairing predicted by 
SMMC is an example of the reentrance (or partial or- 
der) phenomenon manifesting itself in successive phase 
transitions [3J, [35J, |3g, l37|, [38J. One of the signatures 
of the reentrance of the partial order is the anomalous 
specific-heat behavior [35j. Figure [4] shows C v at differ- 
ent rotational frequencies. At small values of to, a typi- 
cal pattern of C v is seen that is characterized by a local 
maximum around fcT=0.6 MeV associated with breaking 
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FIG. 4: (Color online) Specific heat calculated in SMMC for 
72 Ge as a function of temperature at four values of lu (in 
MeV). At w=0.5MeV, C v sharply decreases locally, then in- 
creases and decreases again; this is a signature of the reen- 
trance behavior. 

of pairs [ll[ . This maximum is shifted to slightly lower 
temperatures at higher rotational frequencies, due to the 
Coriolis-antipairing effect. At the lowest temperatures, 



C v shows a sharp rise, absent at uj=0 [llj. This can be 
attributed to a nonzero alignment (cf. Fig. [3]) that af- 
fects the nuclear moment of inertia. At w=0.5MeV, the 
specific heat exhibits a clear local dip associated with the 
reentrance of pairing. 

In summary, using the SMMC technique, we explored 
properties of 72 Ge as a function of temperature and 
polarization, paying particular attention to the nuclear 
Meissner effect and dynamics of angular momentum gen- 
eration. We find that a significant contribution to the 
average angular momentum of a heated nucleus comes 
from quantum and thermal fluctuations of the angular 
momentum in the direction perpendicular to that of the 
average polarization. Our calculations demonstrate the 
presence of the partial order phenomenon associated with 
the reappearance of pairing at high rotational frequencies 
and intermediate temperatures. 
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